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ABSTRACT 

Properties of dark matter haloes can be probed observationally and numerically, and comparing both approaches provide ways to 
constraint cosmological models. When it comes to the inner part of galaxy cluster scale haloes, interaction between the baryonic and 
the dark matter component is an important issue which is far to be fully understood. With this work, we aim to initiate a program 
coupling observational and numerical studies in order to probe the inner part of galaxy clusters. In this article, we apply strong lensing 
techniques in Abell 1703, a massive X-ray luminous galaxy cluster at z = 0.28. Our analysis is based on imaging data both from space 
and ground in 8 bands, complemented with a spectroscopic survey. Abell 1703 is rather circular from the general shape of its multiply 
imaged systems and is dominated by a giant elliptical cD galaxy in its centre. This cluster exhibits a remarkable bright 'central ring' 
formed by 4 bright images at Zspec = 0.888 only 5 - 13" away from the cD centre. This unique feature offers a rare lensing constraint 
to probe the central mass distribution. The stellar contribution from the cD galaxy (~ 1.25 10'- Mq within 30kpc) is accounted for in 
our parametric mass modelling, and the underlying smooth dark matter component distribution is described using a generalized nfw 
profile parametrized with a central logarithmic slope a. The rms of our mass model in the image plane is equal to 1.4". We find that 
within the range where observational constraints are present (from ~ 20kpc to ~ 210kpc), a is equal to 1.09^jj (3cr confidence 
level). The concentration parameter is equal to C200 ~ 3.5, and the scale radius is constrained to be larger than the region where 
observational constraints are available (r, = 730^,5 kpc). The 2D mass is equal to M (210kpc) = 2.4 lO''' Mq. Although, we cannot 
draw any conclusions on cosmological models at this point since we lack results from realistic numerical simulations containing 
baryons to make a proper comparison. We advocate the need for a large sample of well observed (and well constrained) and simulated 
unimodal relaxed galaxy clusters in order to make reliable comparisons, and potentially provide a test of cosmological models. 

Key words. Gravitational lensing: strong lensing - Galaxies: clusters: individual (Abell 1703) - 



1. Introduction 

Large N-body cosmological simulations have been carried out 
for a decade, with the goal of making statistical predictions on 
dark matter (DM) halo properties. Because of numerical issues, 
most of these large cosmological simulations contain dark matter 
particles only. They all reliably predict that the 3D density pro- 
file Pdm('") should fall as r^^ beyond what is usually called the 
scale radius. Obse rvations have confir med these predictions {e.^. 
Kneib et al. 2003; Pointecouteau et a l. 2005; M andelbaum et al] 
2008). This agreement is likely to be connected with the fact 
that at large radius, the density profile of a galaxy cluster is 
dark matter dominated and the influence of baryons can be ne- 
glected. On smaller scales, if we parametrize the 3D density 
profile of the DM using a cuspy profile pdm °^ dark mat- 



ter only simulations predict a logarithmic slope a ~ 1-1.5 for 
r — > 0. The exact value of the central slope and its universality 
is deb ated dN avarro et aOl997tlM oore et alll998tlGhigna et all 
2000'; Ricotti 2003i:iNavarro et al.ll2004t iGaoet alJ l2008\ From 
the theoretical point of vie w, the logarit hmic s lope of DM halos 
is predicted to be a ~ 0.8 dAustin et a l. 2005; Hansen & Stadell 
2006). Although this debate is of some interest, these dark mat- 
ter only studies and their predictions do not help much to make 
comparison with observations. Indeed, observing the central part 
(i.e. the inner ~ 500 kpc) of a galaxy cluster at any wavelength 
reveals the presence of baryons (in the forms of stars and X-ray 
hot gas). Thus any attempt to compare observations to simula- 
tions in the centre of galaxy clusters has to be made with numeri- 
cal simulations (or calculations) taking into account the baryonic 
component and its associated physics. 
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Efforts are currently developed on the numerical side in order 
to account for the presence of baryons (e.g. the Horizoifl simu- 
lation). In practice, our understanding of the baryonic physics is 
poor, and the exact interplay between dark matter and baryon is 
far from understood. For example, due to the overcooling prob- 
lem, numerical simulations predict blue central brightest clus- 
ter galaxies which are not always observed. Different effects do 
compete when it comes to the central slope of the density pro- 
file: the cooling of gas in the centre of dark matter haloes is 
expected to lead to a more concentrated dar k matter distribu- 
tion (the so-called adiabatic contraction, s ee Blumenthal et al.l 
119861; iGnedin et al.ll2004 iGustafsson et al.ll200 6'). On the other 
hand, dynamical friction heating of massive galaxies against 
the diffuse cluste r dark matter can flatten the slope of t h e DM 
density profile (El-Zant etalj|200lt iNipoti et alj l2003l l2004t 
iMa & Bovlan' Kolchin 20047 and this effect could even domi- 
nate over adiabatic contraction -Zant et al. 1 12004). Note also 
that the properties of the inner part of simulated galaxy clus- 
ters (even in dark matter only simulations) c an depend signifi- 
cantly on initial conditions as demonstrated in lKazantzidis et aU 
( |2004|) . To summarize, no coherent picture has yet emerged from 
N-body simulations when it comes to the shape of the inner den- 
sity profile of structures; this problem is a difficult one and the 
answer is likely not to be unique but may depend on the physical 
properties of the structures and their formation history. 

On the observational side, efforts have been put on 
probing the central slope a of the underlying dark mat- 
ter distribution. These analyses have led to wide-ranging re- 



lArabadiis et al."2002 


' — - — — — ■ — ' 

; Lewis et al. 2003; Zappacosta et al. 200^; 


lensing (Tvson et al. 


1998; Smith et al. 2001; Dahle et al. 20031 



2004, 2008; Bradac et a l..2007l) or dvnamics (Kelson et al..2002 



Biviano & Salucci 200^. This highlight the difficulty of such 
studies and the possible large scatter in the value of a from one 
cluster to another. 

To summarize, we need to probe observationally and numer- 
ically the behaviour of the underlying dark matter distribution 
(i.e. after the baryonic component has been separated from the 
dark matter component) in the central parts of galaxy clusters. 
The main difficulties are: i) Observationally, to be able to disen- 
tangle the baryonic component and the underlying dark matter 
distribution; ii) Numerically, to implement the baryonic physics 
into the simulations; iii) Then to compare both approaches in 
a consistent way. These issues are far beyond the scope of this 
article, but they are likely to provide an interesting test of the 
ACDM scenario in the future. 

In this work, we aim to probe observationally the central 
(i.e. from ~ 5" (~ 20kpc) up to ~ 50" (~ 200 kpc) from 
the centre) density profile of a massive cluster lens. Our main 
goal is to measure the slope of the inner underlying dark mat- 
ter distribution within this radius. In practice, we apply strong 
lensing techniques on galaxy cluster Abell 1703, a massive 
z = 0.28 (Allen et al. 1992) X-r ay cluster with a luminosity L^ 
= 8.66 lO'^'* ergs ' dBohringer et al. 2000). It is very well suited 
for the analysis we want to perform for the following reasons: 

1. It contains a large number of gravitational arcs, providing 
many lensing constraints for the analysis (Fig.[rii. 

2. Although it displays an intriguing filamentary structure 
along the north-south direction, it looks rather circular from 
the geometrical configuration of its multiply imaged sys- 
tems, in particular its giant arc, located at large angular sep- 



aration (~ 35" ~ 147 kpc). It is likely to be a unimodal 
cluster, with a clear dominant elliptical cD galaxy, which 
makes it much easier to interpret the results of the mod- 
elling compared to bi modal clusters such as Abell 168 9 
jLimousin et al.'20 07bl) Abel l 22 1 8 ( Elfasdottir et al.ll2007l) . 
Abell 68 iRichard e Tal]l2007l) or MS 2053.7-0449 (Verdugo 
et al., 2007). Indeed, regular relaxed clusters are rare at such 
redshifts (Smith et al. 2005|). 
3. It presents a remarkable lensing configuration, forming a 
'central ring' composed of four bright images (see Fig.[T]and 
Appendix). Interestingly, these constraints are found very 
close to the cD galaxy, potentially providing robust con- 
straints in the very central part of the cluster. Note that this 
'central ring' is not really central since it does not coincide 
with the peak of the total mass distribution. 

This article is organized as follows: data used in this work are 
presented in Section 2. Section 3 presents the multiply imaged 
systems that constitutes the basis of our analysis. The strong 
lensing analysis is outlined in Section 4, and the results are given 
in Section 5. We discuss our results and conclude in Section 6. 

All our results are scaled to the flat, ACDM cosmology with 
Qm = 0.3, Oa = 0.7 and a Hubble constant Hq = 70 kms"' 
Mpc"'. In such a cosmology, at z = 0.28, 1" corresponds to 
4.244 kpc. All the figures of the cluster are aligned with WCS 
coordinates, i.e. north is up, east is left. The reference centre 
of our analysis is fixed at the cD centre: R. A. = 13: 15:05.276, 
Decl.=+5 1:49:02.85 (J 2000.0). Magnitudes are given in the AB 
system. 
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Fig. 1. Colour image of Abell 1703 from F850W, F625W and F475W observations (| Stott|l2007 f). North is up, east is left. Size of the 
field of view is equal to 77" x 107", corresponding to 326 kpc x 454 kpc. Multiply imaged systems used in the analysis are shown. 
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Table 1. Multiply imaged systems considered in this work. We have found 13 distinct multiply imaged systems. Coordinates are 
given in degrees (J 2000.0). When the photometry is reliable in each band, we report the corresponding photometric redshift esti- 
mates, with error bars quoting the 3cr confidence level. A(zphot) corresponds to the redshift range allowed by the photometric redshift 
estimation and that will be used as a prior in the optimization (i.e. for each system, the redshift will be let free and allowed to vary 
between A(zphot)). For system 1 however we fix the redshift to the measured one, 0.8885. The photometric estimate A(zphot) - [0.725 
- 1 .04] is in agreement with the spectroscopic measurement. Zmod corresponds to the redshift inferred from the optimization pro- 
cedure. We report both the RMS in the source plane and the RMS in the image plane. The mean scatters are given for the whole 
system, not for each individual image composing a system. The total RMS is equal to 0.26" (source plane) and 1 .45" (image plane). 
Zfree correspouds to the redshift inferred from the optimization when all redshift but system 1 are assigned a flat prior between 0.28 
and 6 (see Section 5.3). 
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Table 2. Different observations of Abell 1703 used in this work: 
instrument, filter and exposure time in seconds. 



2. Data 

Abell 1703 has been observed from space by the Hubble 
Space Telescope (using ACS and NICMOS); from ground 
with the Subaru telescope (using th e Multi-Object Infra Red 
Camera and Spectrograph, MOIRCS. Ilchikawa et'ani2006l) and 
with the Keck telescop e (using the Lo w Resolution Imaging 
Spectroscograph, LRIS. lOke et all 1 19951) . Table |2] sums up the 
different dataset. 

2.1. Imaging Data 

The multiwavelength ACS data have been used to construct 
colour images of Abell 1703 in order to identify the mutiply im- 
aged systems. In addition to the 6 ACS bands, we benefited from 
NICMOS and Subaru data in order to construct spectral energy 
distributions (SED) and estimate photometric redshifts for the 
multiple images as well as the stellar mass of the cD galaxy. 

HST/ACS: We use the HST ACS data in six bands. Abell 1703 
has been observed on May 2005 as part of the ACS team guaran- 
teed observing time (RI. Holland Ford, proposal 10325). These 
images have been reduced using the multidrizzle software to 
remove cosmic rays, bad pixels, combine the dithered frames 
and correct for geometric distortions. The output pixel scale was 
fixed at 0.04"and we used a pixfrac parameter value of 0.8 for 
reducing the area of the input pixels. 

HST/NICMOS: Abell 1703 has been observed with NICMOS 
on October 10, 2006 (RI. Holland Ford, proposal 10996), using 
the FllOW filter. NICMOS images have been reduced follow- 
ing the NICMOS data reduction handboolfl and including spe- 
cific improvements to remove cosmic rays, quadrant-to-quadrant 
variations, bias and flat residuals. More details on these improve- 
ments are given in Richard et al.l (i2008.) . 

Subaru MOIRCS: MOIRCS observations of Abell 1703 were 
obtained in May 2007 (PI: E. Egami). We took five series of 
dithered exposures, with individual exposure times of 20, 25 and 
30 seconds and 5 coadds per frame, producing a total exposure 
time of 18 875 seconds in the central region. These observations 



were obtained under excellent seeing conditions (0.3-0.4"). We 
used the MCSRED packag^ to perform flat-fielding, sky sub- 
traction, distortion correction and mosaicking of individual im- 
ages. The photometric calibration was derived with the observed 
magnitudes of 15 stars from the 2MASS catalog. 

2.2. Spectroscopic Data 

We used LRIS on Keck I in an attempt to measure a redshift for 
the brightest component of system 1 (image 1.1). Four exposures 
of 900 seconds were taken on Jan 29th 2008, under photometric 
conditions but a poor seeing (1.4"), using the 900 lines. mm"' 
grating blazed at 6320 A in the red channel of the instrument. 
This setup covers the wavelength range 5600-7200 A at a resolu- 
tion of 2.77 A and a dispersion of 0.83 A per pixel. The spectrum 
has been reduced with standard IRAF procedures for bias cor- 
rection, flat-fielding, sky subtraction and distortion correction, 
in that order We used the numerous sky lines in this region for 
the wavelength calibration and observation of the standard star 
Feige 92 for the flux calibration. 

We detected a bright doublet of emission lines, centered 
at 7036 and 7041 A respectively, which we interpret as 
[0//]3726,3729 at z = 0.8885 ± 0.0002 without any ambiguity, the 
doublet being easily separated by 5A at this resolution (Fig.|2]i. 
The spectroscopic redshift is in agreement with our photometric 
redshift estimate of Zphot = 0.965!|];^^f) (Table©. 

2.3. Photometric Analysis 

We created photometric catal ogs combining the multi colour im- 
ages by running SExtractor (iBertin & Arnoutslll996h in double 
image mode. A specific detection image was created by combin- 
ing all the ACS data after scaling them to the same background 
noise level. We used the cD-subtracted images (see Section l2!6l ) 
for ACS and MOIRCS to prevent any strong photometric con- 
tamination in the central regions, in particular for the measure- 
ment of the lensed ring-shape system (system 1). Photometry 
was optimized to the small size of the lensed background galax- 
ies, by measuring the flux in a 1.0"diameter aperture. This gives 
accurate colours in the optical bands that used the same instru- 
ment. We estimated aperture corrections in the near-infrared by 
measuring the photometry of 10 isolated bright point sources in 
the NICMOS and MOIRCS images, and corrected our photo- 
metric catalogs for this difference. 

We increased the photometric error bars, usually underes- 
timated by SExtractor, to take into account the effects of driz- 
zling in the reductio n of ACS and NI CMOS images, following 
file computations by ICasertano et al.l (2000). For the MOIRCS 
images, we measured the pixel-to-pixel background noise from 
blank regions of sky selected in the original images, and scaled 
it to the aperture size used in the photometry. 

2.4. Cluster Member Identification 

To extract cluster galaxies, we plot the characteristic cluster 
red sequences (F775W-H) and (F625W-H) in two colour- 
magnitude diagrams and select the objects lying on both red- 
sequences as cluster galaxies. This yields 345 early-type cluster 
galaxies down to F775W = 24. For the purpose of the modelling, 
we will consider only the cluster members whose magnitude is 
brighter than 21 in the F775W band (in order to save comput- 



^ http://www.stsci.edu/hst/nicmos/documents/handbooks 



^ http://www.naoj.org/stafF/ichi/MCSRED/mcsred.html 
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Fig. 2. Spectroscopic observation of image 1.1. Left: Position of the slit. Right: ID spectrum. The [C//]3726,3729 doublet is easily 
identified. 



ing time, see lElfasdottir et ani2007h and which are located close 
to some multiply imaged systems. We consider fainter galaxies 
only if they are located close to some multiple images sin ce they 
can lo cally perturb the lensing configuration (Meneghetti et alj 
l2007al) . This yields 45 galaxy scale perturbers. 



2.5. Photometric redsiiifts 



We ran the photometric redshift code HyperZ jBolzonella et alj 
^000) on the multi-band photometric catalogs, in order to get 
a redshift estimate for all the multiple images identified in 
the field. This program performs a minimization procedure 
between the spectral energy distribution of eac h object and 
a library of s pectral templates, either empirical (IKinney et al.l 
19961; Coleman et al. 1980) or from the evolutionary models by 



Bruzual & Chariot! (l2003l) . Note that the NICMOS data we have 



correspond to a mosaic of four pointings, none of them being 
centred on the cD galaxy. It results that the cD galaxy is not 
fully sampled by the NICMOS data. Therefore, it cannot be sub- 
tracted in the NICMOS image, thus we did not take into account 
the FllOW filter measurement for all its neighbouring objects, 
or for multiple images located at the edges of the NICMOS field 
of view. Absolute photometric calibration between the different 
bands is usually accurate down to 0.05 magnitudes, we used this 
value as a minimal value in the code (this value will be used 
also for fitting the SED of the cD galaxy in Section 2.6). We 
scanned the following range of parameters: 0.0 < z < 7.0 for 
the redshift, 0.0 < Ay < 1 .2 for the reddenin g, applied on the 
template spectra with the ICalzetti et al.l ([2000) law observed in 
sta rburst galaxies . Optical depth in the Lyman-a forest followed 
the lMadaul([T995l) prescription. 

The results obtained for each multiple image are reported in 
Table [T] along with the 3cr error bar estimate given by HyperZ 
from the redshift probability distribution. 



2.6. cD galaxy 

Abell 1703 exhibits a dominant central giant elliptical galaxy: 
its stellar contribution to the mass budget in the central part is to 
be taken properly into account. We worked out some of the cD 
galaxy properties from the broad band photometry. In particular, 
we are interested in computing its luminosities in the different 
filters in order to estimate its stellar mass. Note again that we 
have not been able to study the cD galaxy in the NICMOS band 




Distance from center (orcseconds) 

Fig. 3. Integrated luminosity profile of the cD galaxy, in the B 
band rest frame, for the fit with redenning (solid) and without 
redenning (dashed). The corresponding mass profile used in the 
strong lensing analysis is shown as dot-dashed line. Both profiles 
have the same behaviour for R > ~ 2". 



since the NICMOS data correspond to a mosaic of four point- 
ings, none of them being centred on the cD galaxy. 

Subtraction of the cD galaxy: We fitted and subtracted from 
each image a model representation of the surface brightness dis- 
tribution using the IRAF task ellipse. Both the position angle 
and ellipticity were allowed to vary as a function of the semima- 
jor axis in the fitted elliptical isophotes, as well as the isophote 
centroid in the central part. This procedure was found to give 
satisfactory residuals at the centre (see Appendix). This proce- 
dure also allows us to determine accurate photometry for the cD 
galaxy, that we report in Table [3] The integrated luminosity pro- 
file in the rest frame B band is shown in Fig. [3] together with a 
mass profile of the cD galaxy as included in the modelling (see 
Section|4!2l). 

Stellar mass: From the broad band photometry, we estimate 
the stellar mass of the cD galaxy by fitting its spectral energy 
distribution using the HyperZ software. In order to perform this 
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Lambda(A) 



Fig. 4. Results of the SED fitting procedure. Luminosities are 
computed in an aperture of 7". The red soHd Hne corresponds 
to the fit wit redenning, and the blue dotted-dashed line corre- 
sponds to the fit without redenning. In each case, the fit is good 
(X^ ^ J < 1) and leads to the same estimation of the stellar mass. 



Filter 


L (no Av) 


L(Av) 


M./L (no Av) 


M./L (Av) 


B 


1.1 10" 


1.610" 


11.2 


8.1 


555 


1.5 10" 


2.210" 


7.8 


5.9 


606 


1.710" 


2.5 10" 


6.8 


5.2 


775 


1.510" 


2.210" 


7.7 


6.0 


H 


7.210" 


9.5 10" 


1.6 


1.4 



Table 3. Luminosities and stellar mass to light ratio computed 
in different rest frame filters. The B filter corresponds to the stan- 
dard Bessell B filter, and values are given in this filter in order to 
allow easier comparisons with other galaxies studied by differ- 
ent authors. All values are given in solar units. We report values 
corresponding to two different fits, with and without reddening 
respectively. 



estimation, we choose an aperture equal to 7". This choice is 
motivated by the goodness of the fit in this aperture (x^ ^ ^ < 1). 
Moreover, we find the stellar mass to light ratios (M»/L) to be 
constant within this aperture. We use 8 evo lutionary synthetic 
SEDs computed w ith the las t versi on of the iBruzual & Charlon 
(l2003h code, with Ch abrieJ (12003') IMF and solar metallicity. 
Given the data available, the photometric SED is equally well 
fit with either a very old elliptical template (a Single Stellar 
Population, SSP, aged 10.2 Gyr), or a slightly younger popu- 
latio n (SSP aged 6-7 G yr) including some reddening (Ay =0.4, 
with lCalzetti et alj (1200 0) extinction law). The degeneracy in the 
(Ay, age) plane yields slightly different results in the best fit tem- 
plate depending on the aperture radius. However, both types of 
models (i.e. with and without reddening) are equally likely, as 
shown in Fig. |4] for the adopted 7" aperture. Although the total 
luminosity and M^/L ratio depend on the reddening correction, 
the stellar mass within a given aperture is found to be insensi- 
tive to the best fit template. The resulting stellar mas s within a 
7" aperture is M. = 1.25 ±0.3 lO'^M©. According to'B ell et alj 
(12007), the errors coming from setting the overall mass scale and 
its evolution are usually 25%, much larger than the error on the 
photometry. We adopt 25% accuracy on the stellar mass estima- 



tion. Given the luminosity of this galaxy, we derive a a stellar 
mass to light ratio in the rest frame B band equal to ~ 8 - 11 (so- 
lar units), depending on the redenning. We report all these values 
in Table [3] All luminosities are given in solar units calculated in 
the rest frame filters. 

This is comparable (though a bit higher, but note that this 
galaxy has a very massive stellar population) to typical val- 
ues of stellar mass to light ratio for giant elliptical galaxies 
(iGerhard et al.ll200l . 



3. Multiply Imaged Systems 

Fig. [T] shows a colour image of Abell 1703 where we label the 
multiply imaged systems. We report their positions and photo- 
metric redshifts in Table [T] The identification is a difficult step 
which is done in an iterative fashion: we begin to build a model 
using the most obvious lensed features. Then this model is used 
to test and predict possible multiply imaged systems. In total, we 
use 13 multiply imaged systems in this analysis. It is certain that 
more systems have to be found within the ACS field that presents 
many likely blue lensed features. 

Here we give some notes on the different systems, and 
present colour images for each system in Appendix. 

- System 1, the 'central ring': it is composed by 5 images. The 
four main bright images located close to the cD galaxy dis- 
play a rare "hyperbolic umbilic" lensing configuration (see 
Fig. [1] and Appendix). A demagnified counter image is lo- 
cated in the east. Image 1.3 is located at 5.5" = 23kpc from 
the centre of the cluster, whereas image 1.1 is located at 
13.5" = 56kpc from the cluster centre. This system consti- 
tutes the innermost observational lensing constraint available 
in this analysis. This is the only system for which we have a 
spectroscopic redshift, equal to 0.8885. 

- System 2: This system is constituted by two images forming 
a straight gravitational arc. Two counter images are predicted 
to be more than two magnitudes fainter, we have not been 
able to detect any. 

- System 3: System 3 is located in the northern part of the ACS 
field, close to Galaxy 852. It is composed of two merging 
bright images, and an additional fainter counter image a bit 
further west. 

- System 4-5: We propose that a single lensed galaxy can be 
resolved into two parts, each part being considered as a mul- 
tiply imaged system. It constitutes a typical cusp configura- 
tion system, with the images forming only on one side of the 
cluster. 

- System 6: This is another cusp configuration system located 
a bit closer to the centre than system 4-5. 

- Systems 7-8-9: These three cusp configuration systems are 
located in the north, at a radius a bit further from the centre 
than system 4-5. 

- System 10-11: These two systems constitute the giant tan- 
gential arc. They correspond to two different spots we have 
identified on the giant arc that shows many substructures. 

- System 15 and 16: These two systems composed by 4 im- 
ages 'follow' each other (Fig. [TJ: 15./ and 16./ (/=1,2,3,4) 
are found close to each other, presenting an Einstein cross 
configuration. 

- Radial Arc: We report a radial feature coming out from 
the cD galaxy, composed by two spots (Fig. lA.U . Due to 
the presence of the cD, the estimated photometric redshift 
(z ~ 1 .4) is uncertain, and we have not been able to detect 
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possible counterimages. This likely radial arc is not used in 4. Strong Lensing Analysis 

the analysis however. 

4.1. Methodology 

To reconstruct the mass distribution in Abell 1703, we use a 
parametric method as implemented in the publicly available 
LENSTOO 10 software ( Julio et al. 2007). We use the observational 
constraints (positions of the multiply imaged systems) to opti- 
mize the parameters used to describe the mass distribution: this 
is what we refer to as optimization procedure. The strong lensing 
m ethodology used in this analysis has been described in details 
in [Limousin et al.l (i2007b). We refer the interested reader to this 
article for a complete description of our methodology. 

We describe the mass distribution in Abell 1703 by con- 
structing a two components mass model: the contribution from 
the dominant central cD galaxy is fixed by its stellar mass, and 
then the remaining mass is put into an underlying smooth dark 
matt er distribution de scribed using a generalized NFW profile 
(see ISand et al.ll2008l for details on the implementation of the 
generalized NFW profile into the lenstool code). We also take 
into account the perturbations associated with the galaxies. 



4.2. Modelling the cD contribution 

Degeneracies can arise between the two mass components: if too 
much mass is put into the cD galaxy, then this can lead to a shal- 
lower slope of the dark matter halo and vice versa. Therefore, 
special care has to be taken when modelling the cD galaxy, 
and such a modelling must be as much as possible 'observa- 
tionally motivated' . We use a dual Pseudo Isoth ermal Elliptical 
Mass Distribution (dPIE, see lEUasdottir et al.l20 07) with no core 
radius to describe the cD stellar mass contribution. This pro- 
file is formally the same as the Pseudo Isoth ermal Elliptical 
Mass Distribution (PIEMD) pro file described in 'Limous in et al.l 
(l2005h . However, as explained in lElfasdottir et a l. (2007) it is not 
the same as the PIEMD originally defined bv^Kassio la & Kovnerl 
(il993). Therefore we have adopted the new name dPIE to avoid 
confusion. The position of this mass clump is fixed at (0, 0); the 
ellipticity and position angle are set to be the one we measured 
from the light distribution. Given this parametrization, the mass 
of the cD scales as Mcd oc o-y x rcut, where the scale radius rcut 
is almost equal t o the half mass radius, and cTq is a fiducial ve- 
locity dispersion (lEliasdottir et al.l2007l) . The choice of the scale 
radius will influence the shape of the mass profile. The smaller 
the scale radius, the steeper the mass profile. We choose a scale 
radius of 25 kpc so that the mass profile is as close as possible 
to the integrated luminosity profile. As we can see on Fig. [3] 
the mass profile and the luminosity profile have the same be- 
haviour for R > ~ 2". Note that 25 kpc is close to the value used 
by Sandetalj ( l2008h to describe the cD galaxies in MS 2137 
(22 kpc) and Abell 383 (26 kpc), where the same mass profile has 
been used. The only free parameter describing this cD galaxy is 
then ctq. This parameter is allowed to vary in a range that is set 
by the choice of the scale radius and the estimation of the stellar 
mass. 



http://www.oamp.fr/cosmology/lenstool/ 
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4.3. Modelling the underlying dark matter distribution 

We assume that the dark matter component can be described us- 
ing a generalized NFW model. Its 3D mass density profile is 
given by: 



PcSc 



(r/r,)-(l + (r/r,)) 



3-Q 



(1) 



with r, the scale radius (note that the scale radius for an NFW 
profile does not h ave the same meaning as the scale radius for 
a dPIE profile, see lLimousin et al] (l2005l) \ For a galaxy cluster, 
this parameter is supposed to be larger than ISOkpc (Section 
4.5), so in the range in radius we probe in this work, we have 
r< and thus the density profile can be approximated by p(r) ~ 

Note that we have subtracted only the stellar contribution of 
the cD galaxy in our modelling. This is consistent with the gen- 
eral picture that the dark ma tter halo of the cluster is also the one 
of the cD galaxy (see, e.g. iMiralda-Escudd [l 9951) . We cannot 
distinguish both haloes. 



4.4. Galaxy scale perturbers 

On top of these two components, we include the brightest clus- 
ter members in the optimization (Section 2.3). This galaxy scale 
component is incorporated into the modelling using empirical 
scaling relations that relate their dynamical parameters (cen- 
tral velocity dispersion and scale radius) to their luminosity, 
whereas their geometrical parameters (centre, ellipticity, posi- 
tion angle) are set to the one measured from the light distribution 
(see Limousin et al. 2007b, for details). This galaxy scale com- 
ponent is thus parametrized by only two free parameters, and at 
the end of the optimization procedure, we get constraints on the 
parameters for a galaxy of a given (arbitrary) luminosity which 
corresponds to an observed magnitude mpiisw - 18.3. 

One single galaxy (labeled 852 in our catalogue) is opti- 
mized individually, in the sense that some of the parameters de- 
scribing this galaxy are allowed to vary instead of being fixed 
by its luminosity. This was necessary to reproduce better the 
geometrical configuration of some images falling close to this 
galaxy (systems 3, 7, 8 and 9). In fact, a visual inspection at the 
image of the cluster shows that this galaxy is very bright and 
extended (Fig.[T]). This galaxy is not representative of the clus- 
ter galaxy population, thus the adopted scaling laws might not 
apply to this object. Moreover, as shown in Appendix, a likely 
lensed blue feature is coming out from this galaxy, suggesting a 
massive substructure. 



4.5. Limits on the parameters 

Each parameter is allowed to vary between some limits (pri- 
ors). The position of the DM clump was allowed to vary be- 
tween + 25" along the X and Y directions. Its ellipticity was 
forced to be smaller than 0.5, since beyond that, the elliptic- 
ity as defined in this work for an NFW profile is no longer 
valid dGolse & Kneibll2002h . Its slope a was allowed to vary 
between 0.2 and 2.0 and its scale radius between 150 and 
750 kpc dTasitsiomi et alJl2004tlDoiag et al.ll2004 . The concen- 
tration parameter C200 was allowed to vary between 2.5 and 9, 
which is large enought to include the expec tations from the most 
recent results from N-body simulations (Neto et aLll2007l) . The 
mass of the cD component was forced to be within the range 
allowed by the stellar mass estimate. Concerning the galaxy 
scale component, we allowed the velocity dispersion to vary 



between 150 and 250 kms"', and the scale radius was forced 
to be smaller than 70 kpc, since we have evidences both from 
obse rvations (N atarajan et al. 1998; Geiger & Schneider 199^ 
Nati u-aian et alj l2002albl: iLimousin et alJj 2007a: Halkola et ap 
2007^ and from numerical simulations (ILimousin et al.,,2007cl) 
that dark matter haloes of cluster galaxies are compact due to 
tidal stripping. 



5. Results 

5.1. Mass Distribution from Strong Lensing 

Results of the optimizatiorQ are given in Table 4. The images 
are well reproduced by our mass model, with an image plane 
RMS equal to 1.4" (0.2" in the source plane). RMS for individ- 
ual systems are listed in Table[T] We derive from our model a 2D 
projected mass within 50" equal to M(50") = 2.4 10'"* Mq. 

Dark Matter Component: The logai-ithmic slope of the 3D dark 
matter distribution is found to be equal to 1.09^Q*j'j (3cr con- 
fidence level). We show in Fig. |7] degeneracy plots between a 
and: the mass of the cD galaxy, the concentration parameter 
C200, and the scale radius r^. We find C200 ~ [3.0 - 4.2] and 
M200 ~ 1 .8 10''' Mq. However, we caution that these values rely 
on pure extrapolation from the strong lensing fit and that these 
quantities should be probed using weak lensing and/or X-ray 
data. Indeed, the NFW scale radius is found to be larger than 
the radius over which we have observational constraints (i.e. 
~ 54"). The ellipticity of the mass distribution is small, with 
a/b = 1.13. The location of the DM clump coincide with the 
position of the cD galaxy. No second large scale dark matter 
clump was needed by the data. From previous experience mod- 
elling bimodal galaxy clusters such as Abell 1689, Abell 2218, 
Abell 68 and MS 2053.7-0449, we are pretty confident that there 
is no need for a second large scale DM clump. 

Galaxy Scale perturbers: The parameters inferred for the 
galaxy population (for an observed magnitude equal to 18.3 in 
the F775 band) are: cry - 210 + 4kms"', and r^ut -66 + 5 kpc. 
We checked that the degeneracies between these two parameters 
are the expected ones, in the sense that they follow constant mass 
lines (Moc ctq x rcut). Galaxy 852 is found to be quite massive, 
withM(15")~5 10'2Mo. 

5.2. Mass & Light: A Relaxed Unimodal Cluster ? 

We compare the mass and the light distribution in Fig. |6] We 
find they do compare well, suggesting that light traces mass in 
Abell 1703. Also shown is the light distribution from the domi- 
nant cD galaxy, which is found to be consistent with the one of 
the overall mass distribution within a few degrees. Moreover, we 
find the centre of the DM halo to be coincident with the centre 
of the cD galaxy. 

At first approximation, these facts suggest that Abell 1703 
is a relaxed unimodal cluster This hypothesis should be investi- 
gated further, in particular with precise X-ray observations since 



^ A parameter file containing all the following information, and 
which can be used with the publicly available lenstool software, is 
available at http://www.dark-cosmology.dk/archive/ This file can be 
useful for making model based predictions, e.g. counter-images of a 
multiple image candidate, amplification and mass map and location of 
the critical lines at a given redshift, and it will be updated. 
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we expect the X-ray emission to be centred on the cD galaxy, 
and to present a small ellipticity. 

Besides, we can draw a line from the southern part of the 
ACS field up to its northern part that pass through the cD galaxy 
and that connects very bright galaxies, much brighter than the 
overall galaxy population (Fig. [TJ. Interestingly, when looking 
at Abell 1703 on much larger scales from SDSS imaging (which 
size is 581" x 809"), the whole cluster galaxy population looks 
rather homogeneous, with no very luminous galaxies as we can 
observe along this filamentary structure. In the south of the ACS 
field, this structure seem to have an influence on the formation 
of the giant arc, breaking its symmetry (Appendix). The per- 
turbation associated with the northern part of this structure is 
more obvious to detect. Indeed, the formation of systems 3, 7, 
8, and 9 (both their existence and their geometrical configura- 
tion) is connected with this extra mass component. Moreover, 
we can see that a blue lensed feature is coming out from Galaxy 
852 (Appendix). The fact that we had to constrain individually 
Galaxy 852 also points out that some extra mass is needed in this 
region and that this galaxy (and possibly the other bright galaxies 
defined by this filament) is not representative of the overall clus- 
ter population. One tentative explanation could be that we are 
observing a galaxy group infalling in the cluster centre. Though 
this scenario would need a devoted spectroscopic follow up of 
the cluster members to get some insights into the velocity di- 
mension of Abell 1703. 

5.3. Reliability of the constraints on a 

One of the main goals of this work is to measure the slope of the 
underlying dark matter component, parametrized by a. We want 
to stress out again that we have assumed that the underlying mass 
distribution can be described using a generalized NFW profile, 
but this assumption may not be correct. We try some tests in 
order to check the reliability of our measurement. 

cD modelling: As discussed before, a will be sensitive to the 
parametrization of the cD galaxy. During the first steps of this 
work, we noticed that if too much mass is put into the cD com- 
ponent, then a can become smaller than what is found here, 
and vice versa: if not enough mass is put into the cD compo- 
nent, a can become larger However, these scenarios were lead- 
ing to stellar masses that are not compatible with the photometry. 
Remain the choice of the scale radius of the cD galaxy. The value 
used in this analysis is set by the shape of the integrated luminos- 
ity profile. We tried to use smaller (down to lOkpc) and larger 
(up to 45 kpc) values of the scale radius, tuning up or down the 
velocity dispersion to keep the stellar mass of the cD consistent 
with what we have estimated from the broad band photometry. 
We find this does not have strong influence on the estimation of 
a. More precisely, when using a scale radius of lOkpc, we find 
a - 1.10 + 0.04, whereas when using a scale radius of 45 kpc, 
we find 0-= 1.05 + 0.05. 

System 1 : system 1 constitutes the innermost observational 
constraint in our analysis. The configuration of the brightest 4 
images forming the 'central ring' allows us to probe the poten- 
tial at small radii and this system constitutes our most stringent 
constrain: if removing system 1 from the analysis, we found that 
the slope was basically unconstrained, since we could get an 
equal fit with a ~ 0.5 or a ~ 1.5, suggesting that the use of 
a generalized NFW profile was not necessary. We note that the 



RMS of system 1 is slightly larger than the mean total RMS. 
This interesting lensing feature is not perfectly retrieved in our 
analysis. One could try to improve the situation by considering 
the influence of galactic substructures on the formation of this 
ring. However, from the large separation of the counter image 
1.5, it is clear that this central ring is not a galaxy-scale lens. 
System 1 is also the only system for which we have been able to 
measure a spectroscopic redshift. This allows us to check how 
the results presented in this work are dependent of the redshift 
of this system. Instead of fixing its redshift to the measured one, 
z - 0.8885, we allowed it to vary between 0.72 and 1.0 as con- 
strained from the photometry. In this case, the model slightly 
underestimate the redshift and find z ~ 0.78. The slope of the 
underlying dark matter halo shifts to higher values, a ~ 1.3. 
This highlight the importance of having spectroscopic redshifts 
over photometric redshifts when it comes to detailed study of 
mass distributions. 

Photometric redshifts: As we can appreciate in Table 1, most 
of the photometric redshifts are well constrained by the avail- 
able filters. In particular, the redshift probability functions were 
dominated by a well defined peak. In order to investigate the 
importance of using reliable photometric redshifts, we redid the 
analysis as follows: the redshift of system 1 was fixed to the 
spectroscopically measured value, and all other redshifts where 
assigned a flat prior between 0.28 (redshift of the cluster) and 
6. Results of the optimized redshifts (zfree) are given in Table 1. 
We can see that they differ significantly from their photometri- 
cally constrained values, and that they are systematically over- 
estimated. Results on the parameters of the dark matter clump 
are the foflowing: a = 1.22+° °^; C200 = 3.4 + 0.3; r, ~ 520 kpc. 
Therefore, mainly the value oiF the slope is significantly changed 
when letting redshifts being free. Since increasing a will lead to 
decreasing the mass enclosed in a given radius, redshifts are be- 
ing overestimated in order to compensate for. This highlight the 
importance of photometric redshift when spectroscopic redshift 
are not available. 

X-ray gas component: In this analysis, we have not been able 
to subtract the mass contribution from the X-ray gas. Recent 
analysis by Bradac et al. (2007) in RXJ 1347.5-1 145, combining 
lensing and X-ray, has been able to disentangle each component. 
Indeed, they found that the total mass profile (DM+gas-Hstars) is 
slightly higher than the DM mass profile, the ratio of both quan- 
tities at 75 kpc being equal to 85%. This could suggest that one 
can neglect the gas component at first approximation. Of course, 
Abefl 1703 and RXJ 1347.5-1145 look very different thus we 
cannot compare them quantitatively. To be more precise, we 
would need devoted X-ray observations of Abell 1703. 

Limits on the parameters: The Hmits adopted in this work are 
motivated by results from numerical simulations. Here we want 
to investigate further what is the influence of the adopted limits 
on the results, in particular on the slope of the dark matter distri- 
bution. First, we released these limits for C2oo- When allowed to 
vary between 1 and 13 (instead of 2.5 and 9), we find the param- 
eters of the dark matter distribution to be fully consistent with 
the results given in Table 4. 

Another concern is the scale radius r,. We find it to be larger 
than the range within which observational constraints are avail- 
able. Since degeneracies arise between the scale radius and the 
slope of the generalized NFW profile, we redid the analysis by 
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imposing the scale radius to be within different limits in order to 
see the influence on the results. Since we do have observational 
constraints up to 210kpc from the centre of the cluster, we are 
confident that if was below 200 kpc, our analysis would have 
been able to constrain it. Therefore, we consider the following 
ranges: (1) 200-300 kpc; (2) 300-400 kpc; (3) 400-600 kpc and 
(4) 600-800 kpc. We report the results in Table 5: (0) corresponds 
to the mass model presented in Section 5.1 and Table 4. The fol- 
lowing lines correspond to models obtained when assuming a 
different range of limits for the scale radius. For each model (/), 
we report the inferred values for (r,, C200 , «) and compare each 
run with model (0). In order to quantify this comparison, we re- 
port A(log(Ev)) which is the difference between the Bayesian 
Evidence of model (0) and the Bayesian Evidence of the consid- 
ered model. When this quantity is positive, it favours model (0) 
as being more likely. We also report A{x^) which is the differ- 
ence between the of model (0) and the of the considered 
model. When this quantity is negative, it favours model (0) as 
being a better fit to the data. 

We can draw the following conclusions: for models (1), (2) 
and (3), we see that the scale radius is always found at the higher 
end of the allowed limit, suggesting that this parameter could be 
larger than the upper limit assigned. As a result, we find these 
models to be less likely than the model (0), since both their x^ 
and Bayesian Evidence are worse. The results for model (4) are 
essentially the same as for model (0): the parameters, as well as 
the x^ are found to be fully consistent. The Bayesian Evidence 
favours model (4) over model (0), which can be understood by 
the fact that the ratio between the posterior and the prior is 
smaller in the case of model (4). Note that model (3) is also 
consistent with model (0), in the sense that the mass clump pa- 
rameters agree with each other, and that both Evidences and x^ 
are comparable. 

These results suggest that, even though the preferred value 
for rs is found larger than the range over which observational 
constraints are found, the multiple images actually are sensitive 
to the value of r,. This can be understood as follows: the gen- 
eralized NEW profile is not a power law model, in the sense 
that the 3D density p{r) is changing at each radius r. This mass 
profile will be close to isothermal for r ~ r^, and the radius 
at which the profile becomes isothermal could be felt by ob- 
servational constraints located at r < r,. To check this sce- 
nario, we compute the 2D aperture masses inferred from each 
model and compare them. If there is a significant mass differ- 
ence between each model at the radius where observational con- 
straints are present, then we could understand why the observa- 
tional constraints are able to discriminate between each model. 
Comparison between masses, expressed as a percentage (calcu- 
lated as (model(0)-model(l))/model(0)), is shown on Fig.|5] We 
see that mass differences can reach up to 3% at the radius where 
the outermost observational constraint is found (~ 210kpc). The 
question is to know whether we are sensitive to such a mass dif- 
ference. In other words: is the accuracy on our mass measure- 
ment below this mass differences ? From the mcmc realizations, 
we estimate the accuracy on our mass measurement and express 
it as a percentage. This accuracy depends on the distance from 
the cluster centre, and is plotted on Fig.|5] If mass differences are 
below the accuracy on our mass measurement up to ~ 125 kpc, 
we see that for R> 150 kpc, mass differences between model (0) 
and models (1) and (2) are above the accuracy on the mass mea- 
surement, which means that we are able to discriminate between 
these different mass models. This shows that the multiple images 
located further away than 150 kpc from the centre are sensitive 
to the value of r,. 
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Table 5. Results of the analysis when using different limits 
for the scale radius r,. Error bars correspond to Icr confidence 
level as inferred from the mcmc optimization. When the posterior 
probability distribution is not Gaussian, we report the mode and 
asymmetric error bars. (0) corresponds to the results of the mass 
model presented in Table 4. The following lines corresponds to 
models obtained when assuming a different prior for the scale 
radius. For each model (/), we report the results on the parame- 
ters of the dark matter clump. To compare the different lines, we 
report A(log(Ev)) and N(X^). 
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Fig. 5. Mass differences between model (0) and model (1) 
(solid), model (2) (dashed) and model (3) (dotted). The dot- 
dashed line correspond to the accuracy on the mass measure- 
ment, expressed as a percentage. If mass differences are below 
the accuracy up to ~ 125 kpc, we see that for R>150kpc, mass 
differences between model (0) and models (1) and (2) are above 
the accuracy on the mass measurement, which means that we are 
able to discriminate between these different mass models. 

Ultimately, the scale radius of Abell 1703 should be con- 
strained by a carefuU strong and weak lensing and/or X-ray anal- 
ysis. 

6. Discussion 

The main result of the presented work is to measure a cen- 
tral slope equ a l to ~ -1.1. This is close to the predictions from 
iNavarro et all (Il997l) . However, we want to stress again that 
comparing these values is not relevant since dark matter only 
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Clump S{x) S{y) e d r (kpc) a ctoo cto (km s ' ) 



NFW 
cD 


-1.0±0.1 
[0.0] 


0.5±0.1 
[0.0] 


0.12+0.004 
[0.22] 


64.5+0.3 
[52] 


797 n+U.«4 i nn+O.OlS t i +().J4 

— /O. / / — U.UJo — u.u / 

[30] 


299.4+035 


Galax 852 


[19.0] 


[54.0] 


[0.11] 


[65.5] 


97.1±1.2 


319.5±4.1 


L* elliptical galaxy 










65 9+" '* 


207.8+||o 



Table 4. IVIass model parameters. Coordinates are given in arcseconds with respect to the cD Galaxy. The ellipticity e is the one 
of the mass distribution, expressed as - b^/a^ + b^. Error bars correspond to Icr confidence level as inferred from the mcmc 
optimization. When the posterior probability distribution is not Gaussian, we report the mode and asymmetric error bars. Values 
into brackets are not optimized. Note that the meaning of the scale radius reported here is different for the generalized NFW clump 
and the other mass clumps, described using a dPIE profile with no core. 




Fig. 7. Degeneracy plots between a and (from left to right): cD stellar mass, cioo and r,. a corresponds to the logarithmic slope of 
the underlying dark matter mass distribution parametrized by a generalized NFW profile. C200 and correspond to the concentration 
parameter and the scale radius respectively for this mass distribution. Note that the scale radius is found larger than the range within 
which we have observational constraints. 



simulations, by definition, do not take into account the likely 
influence from the baryonic component on the shape of the un- 
derlying dark matter, thus their predictions cannot be reliably 
compared to what is inferred observationally. 

Similar analyses (i.e. lensing analyses aiming to probe 
the central mass de nsity distrib u tion in gala x y clu sters) have 
been carried out b y [Sandet al.l (l2002i |2004 12008^ (see also 
iGavazzi et al]|200l lOavazzi 200l on MS 2137-23). Studies by 
Sand et al. used the measured velocity dispersion pr ofile of the 



cD galaxy as an extra constrain. The first w ork by ISand et al 



(2004) assumed a circular cluster. As shown bv lMeneghetti et al 
(l2007bi) . this assumption was likely to bias the results to- 
wards shal l ower v alues of the central mass density slope. Then 
ISand et all (l2008h redid their analysis using a full 2D lensing 
analysis, taking into account the presence of substructures and 
allowing the clusters for non circularity. They confirmed their 
earlier claims, in particular, they do find evidence for a central 
mass density slope to be smaller than 1. The results presented 
in this work points out towards higher values for the slope of 
the central mass distribution compared to the studies by Sand 
et al. It is worth mentioning that our analysis is close to the 
ones by Sand et al., since we have both used strong lensing tech- 
niques, and moreover used almost the same lenstool software 
(by the time of the studies by Sand et al., th e Bayesian Monte 
Carlo Markov Chain sampler as described by I Julio et al.l (l2007h 



was not available, and they used a parabolic optimization in- 
stead, but this cannot explain the differences). This wide range 
of slopes found from one cluster to another may point out to an 
intrinsic large scatter on this parameter, which may depend on 
the merger history of the cluster. This possible scatter should be 
probed numerically. 

Prospects on Abell 1703: We can summarize some possible 
extensions of the presented work: 

1 . We have left some blue features that are likely to be lensed. 
Therefore, the multiply imaged systems identification is still 
to be improved. In particular, it would be useful to un- 
derstand the radial feature that is coming out from the cD 
galaxy. 

2. Parametric strong lensing is very sensitive to misidentifi- 
cations, thus we are pretty confident of our identifications. 
However, a measure of their redshifts would be very valu- 
able. 

3. Probing the cluster potential on larger scales: our analysis 
allowed the scale radius to vary within a large prior sug- 
gested by numerical simulations. We have investigated how 
the choice of this prior can influence the results on a, finding 
substantial differences between different priors. Measuring 
reliably the scale radius is thus of first importance and should 
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Fig. 6. Mass contours overlaid on the ACS F850W frame. Also 
shown is the light distribution from the cD galaxy (dashed white 
contours), which is found to be consistent in orientation with the 
one of the mass distribution within a few degrees. Note how the 
light and the mass follow the filamentary structure. Size of panel 
is 153" X 153". 



be done combining strong and weak lensing and/or X-ray 
measurements. 

4. Measuring the velocity dispersion of the stars of the cD 
galaxy in order to add dynamical constraints as in the studies 
by Sand et al. 

5. The more relaxed Abell 1703, the stronger the conclusions of 
the work presented here, and also the easier the comparison 
with N-body simulations. We have argued above that we find 
evidence for a relaxed cluster It would be important to check 
further the dynamical state of Abell 1703. In particular, high 
resolution X-ray observations of Abell 1703 could help to 
answer this question. Moreover, X-ray observations would 
also make possible to su btract the gas contr ibution from the 
mass budget, following Bradac et al. I (l2007l) . 

6. Moreover, the spectroscopic study of the motions of cluster 
galaxies provides clues to the dynamical state of the cluster. 
Such data give an insight into the velocity dimension, and 
can reveal if the cluster is undergoing a merger along the 
line of sight or if it is already well relaxed. Of particular in- 
terest will be to probe the velocity difference between the cD 
galaxy and the brightest cluster galaxies defining this fila- 
mentary structure, in order to investigate if we are observing 
an infalling galaxy group. 

7. Related to this group issue is the study of Galaxy 852 based 
on the detection of a nearby arc and the surrounding multiply 
imaged systems. 

Prospects on the Central Mass Distribution of Galaxy 
Clusters: Obviously, we suffer from small number statistics and 
we need to extend this kind of analyses to other clusters. In try- 
ing to gather a sample of clusters for which the presented analy- 
sis could be applied, one should focuss on relaxed unimodal cD 
dominated clusters which present observational constraints as 
close as possible to the centre of the cluster. In that respect, one 



should look for radial arcs. Note, however, that radial arcs might 
preferentially form in clusters with shallow profile. Ideally, se- 
lection should be independent of arc appearances, but the clus- 
ters with radial arcs will indeed give good constraints. 

On the numerical side, we need to study a sample of many 
galaxy clusters containing baryons, and testing different pre- 
scriptions for the baryonic implementation. We have initiated 
such a numerical study on two galaxy clusters, and results will 
be presented in a forthcoming publication. 

Conducting in parallel an observational and a numerical pro- 
gram is a worthy goal: it is interesting by itself to study what is 
going on in the central part of the most massive virialized struc- 
tures since it can provide insights on the interactions between 
baryons and dark matter particles; moreover, it can potentially 
provide an interesting probe of cosmological models. 
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Fig. A.l. System 1, the central ring composed by 4 bright im- 
ages. For its counter image, see Fig. IA.8I below. Top panel: 
colour image from F850W, F625W and F465W observations. 
Size of panel is 10" x 10". Botom panel: F775W image where 
the light from the cD galaxy has been subtracted. We plot the 
tangential and critical lines at z = 0.88 (redshift of system 1). 
We also report a radial feature that is coming out from the cD 
galaxy (see Fig. 1), that is not used in the analysis since we 
have not been able to detect the counter images. Size of panel 
is 16" X 16". 



Fig. A.2. System 2, constituted by two magnified merging im- 
ages. Two counter images are predicted to be more than two 
magnitudes fainter, we have not been able to detect any. Size 
of panel is 6.3" x 6.3". 



3.1 














3.2 



Fig. A.3. System 3, constituted by two magnified merging im- 
ages and a demagnified one a bit further west. Size of panel is 
4" x 4". 



Appendix A: IVIultiply Imaged systems 

In this Appendix we show colour pictures of the multiply imaged 
systems used in this work. The size of each panel is indicated in 
the caption. All images are aligned with the WCS coordinates, 
i.e. north is Up, east is Left. 
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Fig.A.7. Systems 10 and 11 corresponds to two substructures 
identified on the giant arc. Note than more can be defined along 
this giant arc. Top panel: images 10.1, 11.1, 11.2 and 10.2. Size 
of panel is 18" x 12". Bottom panel: counter image of both 
systems, also presenting substructures. Size of panel is 6.6" x 
6.6". 
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Fig. A.8. Systems 15 and 16. Images belonging to system 15 are systematically brighter than images belonging to system 16. First 
panel, north to south: image 1.5, counter image of the central ring (Fig. lA.ll l; image 15.1 and image 16.1. Second panel: images 
15.2 and 16.2. Third panel: images 15.3 and 16.3. Last panel: images 15.4 and 16.4. Size of each panel is 8" x 8". 
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Fig. A.9. Galaxy 852, located in the northern part of the ACS 
field, at a = 198.76342;^= 51.832523. Size of panel is 8.2" x 
8.2". The blue radial feature is easy to detect on this picture. 
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Fig. A.IO. Subaru image of Abell 1703 (H band). Size of panel is 275" x 452". The filamentary structure identified on Fig. 1 can 
be appreciated further on this image. The frame corresponds to the size of Fig. 1. 



